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ABSTRACT 


The TIROS VII Meteorological Satellite, launched 
19 June 1963, was equipped with a medium resolution 
scanning radiometer. One of the spectral regions ob- 
served, the 15 micron region, largely measured the 
emission from the stratosphere. The spatial and tem- 
poral variations of the observed 15 micron emission 
were found to correspond closely to major thermal 
patterns of the stratosphere. The data have been com- 
piled and mapped on a series of ten-day mean charts, 
normally at five-day intervals, covering the 13-month 
period subsequent to launch. The charts comprise 
data from regions of the globe normally inaccessible 
to conventional means and are presented as a guide 
to stratospheric behavior during the period of June 
1963 to July 1964. 
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AN ATLAS OF STRATOSPHERIC MEAN ISOTHERMS 
DERIVED FROM TIROS VII OBSERVATIONS 


1. INTRODUCTION 

The TIROS VII weather satellite, launched 19 June 1963, included a five- 
channel scanning radiometer among its complement of meteorological sensors. 
One of the radiometric channels was sensitive in the 15 micron spectral region 
where atmospheric emission is predominately due to the 15 micron vibration- 
rotation band of carbon dioxide. The primary goal of the 15 micron experiment 
was to examine the suitability of the spectral region for horizon sensing (Bandeen 
et al. , 1963). The advantages inherent in the 15 micron region are: (1) the very 
short unit optical depths of the troposphere minimize the effect of tropospheric 
phenomena, specifically clouds, on the outgoing emission, and (2) carbon dioxide 
is the only major absorber in the 5 to 30 micron region, the region which encom- 
passes the bulk of atmospheric emission, which is present in sensibly constant 
mixing ratios at altitudes of significant emission. 

Implicit in the aforementioned properties of the 15 micron region is the fact 
that variations in observed emission are due primarily to horizontal variations 
in the temperature field at levels of significant emission. It will be shown that 
most of the emission reaching the satellite originates in the stratosphere. Thus 
a mapping of observed emission duplicates, to a large extent, the major thermal 
patterns of the stratosphere. 

While the extension of in situ measurements of upper atmosphere param- 
eters, especially since the International Geophysical Year of 1957, has been 
accelerated, large areas of sparse data exist even in the Northern Hemisphere. 
The Southern Hemisphere, except for localized areas, is a major data gap. 
Therefore , we feel that even relatively unconventional means such as the remote 
measurements by the 15 micron channel of TIROS VII, should be studied to the 
fullest extent to explore the behavior of hitherto forbidden areas. To quote Prof. 
V.P. Starr, ". . . many new conceptions concerning the operation of the at- 
mosphere are to be looked for at the borderline of measurability" (Starr and 
Wallace, 1964). 

The TIROS VII 15 micron data have been compiled into a series of ten-day, 
global mean-temperature charts covering a period of 13 months. The following 
sections will examine in more detail the physical significance of the 15 micron 
sensor and describe the manner in which the charts were constructed. 
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2. 15 MICRON RADIOMETER 


The 15 micron channel of the five-channel scanning radiometer was sensitive 
from 14.8 to 15.5 microns. The field of view of the instrument was five degrees 
between the half-power points. From the mean satellite altitude of 635km, the 
minimum instantaneous area viewed as projected on the earth was about 2400 
sqkm when the optical axis was normal to the surface. 

The optical axis of the radiometer was at a fixed angle of 45 degrees from 
the spin axis of the satellite. The radiometer viewed in two colinear directions, 
180 degrees apart, known as floor side and wall side. At least one side was ex- 
posed to outer space at all times, providing a check point of absolute zero. The 
spinning motion of the satellite at the rate of about lOrpm, combined with the 
orbital motion with the spin axis approximately fixed with respect to inertial 
space, produced a rather complicated pattern of scanning (Allison and Waraecke, 
1964; also Staff Members, 1964, Vol. 1, App. B). The orbital inclination of the 
spacecraft was 58.2 degrees, prohibiting observations poleward of about 60 de- 
grees latitude. 

A detailed description of the components of the radiometer and its method of 
operation is contained in the TIROS VII Radiation Data Catalog and Users* Man- 
ual (Staff Members, 1964). The emphasis here is on the relation of the measure- 
ments to stratospheric temperatures . 

It is possible to decompose the total outgoing radiation into contributions as 
a function of altitude based on the radiative properties of carbon dioxide and the 
spectral response of the radiometer (Nordberg et al . , 1965). Figure 1 shows the 
results for four model atmospheres and two extremes of nadir angle . ^ (h) ex- 
presses the weight given to various altitudes . 

While the weighting curves are rather blunt and broad, it may be seen that 
most of the outgoing energy originates in the region from 10 to 35 km. The ab- 
solute maxima at low nadir angles is at about 23-25 km which corresponds to 
roughly a pressure of 25-35 mb. Figure 1(c) shows that the weighting functions 
shift upward somewhat at larger nadir angles . To reduce the dependence of the 
observations on nadir angle, only measurements made with nadir angles of 40 
degrees or less were used in this study. 

At low nadir angles , more than 96 percent of the total outgoing intensity 
originates at levels above 10km. Cloudiness, therefore, does not greatly affect 
the observed energy. The effect of cloudiness would be the greatest for a trop- 
ical atmosphere at zero degree nadir angle , due to the high tropopause and large 
difference between surface and tropopause temperature. Since the Tropics is a 
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TEMPERATURE (°K) ¥ (W M ‘ 2 KM' 1 ) ¥ (W M 2 KM') 


la) |b) lc) 

Figure 1 . (a) Typical temperature profiles based on proposed supplements to the 
"U.S. Standard Atmosphere 1962" for 60° North summer, 60° North winter (warm 
and cold) and 15° North (b) Weighting functions ^(h), applying to the measured 
outgoing radiance; nadir angle = 0° (c) Weighting functions, ^(h), applying to 
the measured outgoing radiance; nadir angle = 58° 

region of persistent high cloudiness, the effect of clouds was calculated from 
radiative transfer theory to obtain a limiting "worst case" (Kunde, 1966). The 
result was found to be that an overcast at the tropopause would reduce the de- 
rived mean temperature by only three degrees Kelvin as compared with clear 
skies . At higher latitudes , the effect is much less . Due to the very weak hori- 
zontal temperature gradients in the tropical stratosphere, cold pockets associ- 
ated with areas of persistent high cloudiness are frequently found. The lower 
readings may be partially due to actual colder stratospheric temperatures 
caused by cloudy insulation of the stratosphere from tropospheric heat sources. 

The observations at higher latitudes are more physically meaningful in the 
context of stratospheric dynamics partly because the effect of cloudiness is nil 
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but more importantly, because the lapse rates from 10 to 30 km are nearly iso- 
thermal as shown in Figure 1(a) . Weighted mean temperature patterns at high 
latitudes will correspond more nearly to true temperature patterns at a given 
level than in the Tropics where there is a sharp discontinuity at the tropopause. 
An added complication in the Tropics is the double reversal of the horizontal 
poleward temperature gradient from troposphere to lower stratosphere and from 
lower stratosphere to mid-stratosphere . 

The energy received by the radiometer was converted into units of equiva- 
lent blackbody temperature (T BB ) . Equivalent blackbody temperature is defined 
as the temperature of an isothermal blackbody filling the field of view (as in the 
laboratory calibration) which would cause the same response from the radiom- 
eter as does the generally non-Planckian spectral distribution of the radiation 
emerging from the top of the atmosphere in the direction of the satellite. Since 
a vertically-weighted-mean temperature does not easily fit into the framework 
of existing conventional measurements , we feel that the global patterns produced 
from 15 micron data are much more meaningful than the absolute values of equi- 
valent blackbody temperature . The next section describes the manner in which 
the 15 micron data are displayed. 


3. DATA DEFICIENCIES 

While knowledge of stratospheric behavior from existing conventional meas- 
urements was the prime consideration in the format of the data, modifications 
and compromise were necessitated by the following factors: (1) progressive 
degradation of the radiometer, (2) low signal -to-noise ratio of the radiometer, 
(3) errors in geographically locating the data, (4) gaps in global coverage. 

3.1 Degradation 


The radiometer did not have onboard calibration except for a zero degree 
reference when viewing outer space. The performance of the instrument was 
checked by computing mean quasi-global temperatures as a function of time. 
(The quasi-globe is that portion of the earth viewed by TIROS VII, roughly 65 N 
to 65 S latitude . ) It was found that the mean temperatures decreased sharply 
during the first 30 days after launch and decreased at a much slower rate there- 
after. The cause is unknown. In addition, the wall side and floor side recorded 
different readings for the same target viewed a few minutes apart. The differ- 
ence between wall and floor side measurements increased with time and was 
also found to be a function of the temperature of the viewed region. Corrections 
were determined by fitting smooth curves through the global mean wall and floor 
side temperatures and correcting back to temperatures shortly after launch. A 
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detailed description of the degradation is contained in the TIROS VII manual 
(Staff Members, 1964). Briefly each measurement was modified by applying 
a correction, aT, in the form, 


where 


AT f , w (t) = a f w (t) + b f w (t)T f>w (1) 

f = Floor Measurement 

w = Wall Measurement 

a, b = Correction Factors 

t = Time 

T = Uncorrected Observed Temperature 


The correction factors, a and b, were determined from degradation nomo- 
grams in the TIROS VII manual (Staff Members, 1965, Vol. 3). The time scale 
of the degradation was sufficiently slow that the factors could be computed for 
the mid-date of any map and used for all data included in the map, usually that 
observed over a 10-day period. The degradation of the instrument casts some 
doubt on the validity of the absolute value of the measurements but does not 
greatly affect the final temperature patterns. 


3.2 Signal-To-Noise Ratio 

The 15 micron channel measurements were characterized by a random noise 
component with an rms amplitude of about 5 degrees Kelvin. The noise compo- 
nent was largely eliminated by time and spatial averaging which will be described 
later. 


3.3 Geographic Location 


Depending on the angle between the spin axis and local vertical, the radiom- 
eter scans the earth in one of three distinctive patterns (Allison and Warnecke, 
1964). When the spin axis is nearly vertical, the radiometer traces a helical 
pattern on the ground which does not cross the horizon, known as closed mode. 
When the spin axis is nearly horizontal, the radiometer scans the earth alter- 
nately from floor and wall sides, known as the alternating open mode. At inter- 
mediate angles, which constitute the major portion of the total data, the radiom- 
eter describes a series of open arcs, all traced by either the floor or wall side, 
known as the single open mode. 

In the single open mode the horizon is traversed twice during each rotation 
of the satellite , and the known length and time of the scan enable the data to be 
precisely located. In the closed mode, the spin rate of the spacecraft must be 
interpolated from the single open mode data measured just prior and subsequent 
to the closed mode. Since the spin rate varies somewhat due to magnetic torques 
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and other factors, the data are subject to misplacement. In the alternating open 
mode, the computer distinguishes between the alternate floor and wall side meas- 
urements by the respective lengths of the earth-viewed portion of the scans. 
Midway through the alternating mode sequence, the floor and wall scans are ap- 
proximately equal and the computer is apt to mislabel the scans, resulting in a 
gross mislocation of data. 

As mentioned in a previous section, the computer was instructed to reject 
all data with nadir angles greater than 40 degrees to minimize changes in the 
weighting curves. To avoid mislocation of data, a proviso was added that all 
scans with minimum nadir angles occurring within the scan greater than 38 de- 
grees be rejected completely. All ambiguous portions of the alternating mode 
cycle were thus eliminated as well as most of the closed mode . 

3.4 Data Coverage 


The orbital inclination of the satellite, 58.2 degrees, prohibited measure- 
ments poleward of approximately 60 degrees latitude. Limitations in acquiring 
the data from the spacecraft resulted in additional gaps within the confines of 
the quasi-globe. Of the nearly 15 orbits per day experienced by TIROS VII, a 
maximum of nine orbits can be interrogated by the three Command and Data 
Acquisition stations in North America during any given interrogation day. The 
interrogation day includes all orbits interrogated in the series of nine consecu- 
tive orbits which come within range of the acquisition stations during each 24- 
hour period. The interrogation day often occurs on two seperate calendar days. 
The orbits acquired on one interrogation day are compiled on one Final Meteor- 
ological Radiation Tape (FMRT) . 

Two regions of the quasi-globe are never viewed, Central Siberia and south- 
ern South America. On most days, less than the optimum nine orbits are ac- 
quired. In the latter half of the 13 months of mapped 15 micron data, rarely are 
more than four orbits acquired per interrogation day. To fill out the quasi-globe 
and also to overcome random noise, several FMRT's, usually ten, were com- 
bined to produce one map. The rationale for ten-day maps will be covered in the 
next section. 


4. MAPPING PROCEDURES 

While grouping and altering the 15 micron data were necessitated by techni- 
cal deficiencies as outlined in the previous section, ultimately, the character- 
istic spatial and time scales of the major stratospheris phenomena must be con- 
sidered to arrive at a meteorologically useful product. In brief, the stratosphere 
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above 20km is characterized by large scale, slow moving disturbances super- 
imposed on a zonal current. Teweles (1963) has shown that most of the eddy 
activity is explained by zonal wave numbers one through four. Boville (1960), 
using 25 mb temperature data from Resolute, N. W. T. , performed a spectral 
analysis and found that most of the variance is found in the seasonal shift of 
temperature, with a secondary maximum in the 16 to 32 day range associated 
with baroclinic waves . The findings of these authors guided the choice of length 
and time parameters. 

A Mercator map base was selected for the data display as being most suit- 
able for the quasi-globe. The data were gathered on a rectangular grid of equi- 
spaced points, each grid point being assigned the average value of all observa- 
tions falling within its square of influence. The spacing of grid points was 
constant at five degrees of longitude in the zonal direction and variable, ranging 
from five degrees of latitude at the equator to about 2.5 degrees at latitude 60 
degrees, in the meridional direction. The zonal grid spacing permits the reso- 
lution of disturbances up to wave number 36; however, spatial filtering was em- 
ployed to filter such high wave numbers . 

The long time scale of stratospheric phenomena enabled us to compile the 
data into ten-day mean charts at five-day intervals. Actually, the FMRT, rather 
than the calendar, was our basic time unit and the charts comprise 9 to 11 
FMRT's depending on circumstances. During the course of the 13 months of 
observation, there were days or groups of days when no data were recorded. 

The FMRT's included in a chart were then either expanded or contracted to avoid 
overlapping such periods . The days of no data also prohibited the maintenance 
of a constant five day interval between charts throughout the 13 months. The 
term mean chart is somewhat of a misnomer. Certain areas of the quasi-globe 
may have been sampled once, daily, or not at all during the period of a chart. 
Slow moving stratospheric disturbances and strong spatial smoothing are the two 
factors which mitigated the random sampling of data. 

The use of about ten days of data helped to overcome two of the major data 
deficiencies: data gaps and random noise. Throughout much of the 13 months 
of observation, a ten-day period comprises 50 to 70 orbits, and even during 
data-poor periods, 20 to 30 orbits. The average temperature at an individual 
grid point usually comprised hundreds, often thousands of observations. 

After the data had been gathered on the grid array, two factors were dealt 
with before analysis could begin. First, all grid points of the quasi-globe did 
not have an assigned temperature. Any average temperature comprising less 
than 40 observations was rejected because of random noise. Since an electronic 
curve plotter was used for analysis, a complete field of data was required. The 
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grid array was completed by linear interpolation across data gaps in a zonal di- 
rection. This method approximates the technique a human analyst would employ. 

Second, the temperature field contained many short wave perturbations 
which were considered to be noise. Partially, the perturbations were due to the 
quantization of mean temperatures into integer values by the computer. An e- 
lectronic curve plotter cannot selectively ignore such minor variations and 
therefore the temperature was smoothed by a two dimensional spatial filter. 

The filter was designed to leave the longest wavelengths unchanged in ampli- 
tude and phase and eliminate the shortest wavelengths. The filter was construc- 
ted following the method of Wallington (1962) . 

Basically, two different filters each having three weights in the meridional 
direction and nine weights in the zonal direction were passed over the grid array 
and the weighted temperature was assigned to the central value. The same re- 
sult could have been derived in one step using a five-by-seventeen filter but 
would have had to be modified when the operator included external grid points . 
Using two three-by-nine filters in succession resulted in the loss of only the 
northernmost and southernmost line of data. Since the data are cyclic with 
period 2v in the zonal direction, the array was equivalent to a hollow cylinder 
with only the top and bottom lines of data being external points. 

The first filter had a weight array as follows: 
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The second filter was , 
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0.012 
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0.004 
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0.002 
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0.051 

0.058 

0.051 

0.031 
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Both filters are symmetrical in the meridional and zonal directions to eliminate 
any phase shift of the filtered harmonics. By taking the Fourier transform of 
the filter weights, it is possible to compute the response of the filter as a func- 
tion of wavelength (Holloway, 1958). The wavelength response is shown in Fig- 
ure 2 for the equivalent five-by-seventeen filter. Note that for the longest zonal 
harmonics the response is near unity, indicating that they are unaffected by the 
filtering. At a wavelength of three grid intervals, the zonal response is nega- 
tive, indicating a 180 degree phase shift. The response, however, is only 10 
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WAVELENGTH (in grid intervals) 

Figure 2. Amplitude Response of the Spatial Smoothing Operator in the 
Meridional and Zonal Directions as a Function of Wavelength 

percent and disturbances at this wavelength were largely eliminated. Zonal dis- 
turbances at wavelengths of two and four grid intervals were eliminated completely. 

After the smoothing operation, the data were analyzed by an electronic curve 
plotter at intervals of two degrees Kelvin. The resulting charts are presented 
in the next section. 


5. TEN-DAY MEAN CHARTS 

The TIROS VII 15 micron data have been compiled and mapped on a series 
of seventy -three ten-day mean charts, normally at five-day intervals, covering 
the 13-month period from launch (19 June 1963) to 16 July 1964. The charts are 
numbered such that in general a series of successive and contiguous ten-day 
charts (without time overlap) is made up of either the odd-numbered charts or 
the even-numbered charts. However, the phase of the two sets is staggered such 
that the beginning points in time of the maps of say the odd-numbered set occur 
at the midpoints in time of the maps of the even-numbered set or vice versa. 
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Because of occasional interrogation problems resulting in a failure to acquire 
the data successfully, the pattern of time-continuity is broken at several points 
in both sets of charts . The odd-numbered set is the more nearly continuous , 
suffering only one-to-three day breaks surrounding the five entries of "NO DATA 
ACQUIRED" in Table I, plus an eight-day break during the interval 5-12 May 
1964. The even-numbered set suffers five extensive breaks in the intervals iden- 
tified by entries of "NO DATA ACQUIRED" in Table I, plus a seven-day break 
during the interval 10-16 May 1964. Hence, though the final chart is numbered 
78, there are five fewer actual charts represented because of the indicated inter- 
vals where no data were acquired. The degradation-corrected measurements on 
the charts are given in terms of equivalent blackbody temperatures (T BB ) • All 
measurements were made in the nadir angle range 0° - 40° . Other pertinent 
data concerning the charts are documented in Table I. 
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Documentation Data Concerning 15 Micron Charts 
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